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Abstract 

The influence of background states with increasing level of complexity on the strength 
distribution of the isoscalar and isovector giant quadrupole resonance in 40 Ca is studied. It 
is found that the background characteristics, typical for chaotic systems, strongly affects the 
fluctuation properties of the strength distribution. In particular, the small components of the 
wave function obey a scaling law analogous to self-organized systems at the critical state. This 
appears to be consistent with the Porter-Thomas distribution of the transition strength. 
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Nuclear physics has contributed significantly to the recent progress in understanding 
the chaotic aspects of nonlinear dynamical systems, especially in the context of classical 
versus quantum correspondence. The nucleus is particularly well suited for such studies 
due to the intrinsic quantum nature of the nuclear many-body problem on the one hand, 
and the wealth of experimental data on the other. The chaotic nature of the nucleus is by 
now well documented empirically fl[] and seems natural, bearing in mind the many-body 
character of the nucleus and the complicated form of the nucleon-nucleon interaction. 
Still, however, nuclei reveal many regular, collective phenomena and this coexistence 
of chaos and collectivity is a challenging problem for quantitative study 0. Especially 
interesting in this respect are the nuclear giant resonances which carry a large fraction 
of the total transition strength and are located many MeV above the ground state, in 
the energy region which is expected to be dominated by chaotic dynamics. 

The giant resonance, as a short time phenomenon, involves simple configurations 
of one-particle one-hole (lp-lh) type. Chaos may influence the subsequent decay of 
these components which occurs on longer time scales and gradually evolves to more and 
more complex configurations. Eventually, the initial energy deposited in the nucleus is 
redistributed over all available degrees of freedom and the limit of the compound nucleus 
is reached. This is the limit of fully developed chaos and, as a result, quantum stochastic 
methods based on the random matrix theory |3[ or, alternatively, molecular dynamical 
approaches generating chaotic behavior ||, prove appropriate. The process of giant 
resonance formation and its subsequent decay towards the compound nucleus occurs in 
a closed system and the most basic approach is in terms of a single Hamiltonian acting 
in a rich enough Hilbert space such that the relevant degrees of freedom are included. 
This also provides the most natural scheme for the coupling between a collective state 
and the complex background. In quantum mechanical terms one can then speak about 
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the large and the small components of the nuclear wave function 0]. It is the purpose of 
the present paper to work out such a scheme, to identify which ingredients are relevant, 
to study the role of chaos on the giant resonance physical observables and finally, from 
a more general perspective, to contribute to the understanding of the universal aspects 
of the coexistence between collectivity and chaos in small many-body quantum systems. 
We start from the nuclear Hamiltonian which, in second quantized form, reads as 



H = Y^^Ui + 7} V,,jJ,.l<! ;(////'//,. 



,t„t, 

ij,kl 



The first term denotes the mean field which represents its most regular part || while 
the second term is the residual interaction. Clearly a diagonalization of H in the full 
Hilbert space of n-particle n-hole excitations is numerically prohibitive, and may not 
be necessary. Our recent study f?J of level statistics in the subspace generated by two- 
particle two-hole (2p-2h) excitations has shown that it is sufficient to diagonalize the 
Hamiltonian within a truncated subspace of lp-lh and 2p-2h states ||: 

|1) = ala h \0); |2) = a^a^a^a^jO) (2) 

For the present discussion it is more transparent to prediagonalize v in the lp-lh and 
2p-2h subspaces such that 

|1) = £^11); |2> = X;cf|2> (3) 
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Then the Schrodinger equation takes the following form 
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where Ei and E^ denote the energy eigenvalues in the lp-lh and 2p-2h subspaces re- 
spectively and = J2u C\ (l|u|2)C| mediates the coupling between these two spaces. 



The solutions of Eq. (H) yield the transition strength distribution Sf{E) in response to 
an external field F = J2a F 

S F (E) = J2\X[(0\F\1)\ 2 5(E-E N ) (5) 

i,N 

with its energy moments being defined as 

(E n ) = JdES F (E)E n / JdES F {E). (6) 

For the specific case of the quadrupole response in 40 Ca, considered here, we have 
chosen the mean field and residual interaction as in ref. || including two major shells 
above and below the Fermi level. Our study is based on an explicit diagonalization in 
Eq. (4) which involves more than 11,000 states out of which only 26 are lp-lh states. 
Computational restrictions require to limit the number of 2p-2h states. It turns out 
that, including those up to 50 MeV excitation energy, suffices for a realistic description 



of the measured response function | i0|j . This yields altogether 3014 2 + states which 
is numerically manageable. The results displayed in Fig. 1 and Fig. 2 yield a mean 
excitation energy (E) of 30.84 MeV for the isovector and 24.45 MeV for the isoscalar 
transitions, independent of the mixing with 2p-2h states. The latter can be easily 
understood by writing the mean energy as (E) = (0\F^[H, F]\0) and realizing that 
|0) is given by a Slaterdeterminant of occupied single-particle states. For the energy 
dispersion, a = ((E 2 ) - (E) 2 ) 1 / 2 , we obtain 3.41 MeV and 4.28 MeV when only lp-lh 
excitations are considered. Allowing for the coupling to 2p2h-states these values are 
increased to 5.40 MeV and 5.42 MeV, respectively, independent of whether effects of 
the residual interaction in the 2p-2h subspace is included or not. The latter can be 
deduced from the form of |0) and the fact that (E 2 ) = -(0\[H, F^][H, F]\0). Motivated 
by our previous results 0] on the level fluctuations in the prediagonalized 2p-2h space 
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we distinguish three cases: (1) no residual interaction in this space (the corresponding 
strength distribution is shown in Figs, lb and 2b) for which the nearest-neighbor spacing 
distribution is sharply peaked near zero, because of degeneracies, (2) the inclusion of 
particle-particle and hole-hole two-body matrixelements (Figs, lc and 2c) which results 
in a Poissonian distribution, characteristic of the known universality class of generic 
integrable systems |yj (3) the use of the full residual interaction (Figs. Id and 2d) which 
yields the fluctuations of the Gaussian orthogonal ensemble (GOE) |jj characteristic of 
classically chaotic systems. All three cases introduce significant modifications of the 
lp-lh 'doorway' strength distribution (notice the change in magnitude of the transition 
matrixelements) resulting in a gradual reduction of the large components accompanied 
by a simultaneous amplification of the smaller ones when going from case (1) to case 
(3). The isovector response is affected more than the isoscalar one especially when going 
from (2) to (3) (Figs, lc and Id and Figs. 2c and 2d). 

Of course, what is important for the degree of mixing between the spaces spanned 
by |1) and |2) is not only the spectral properties of the eigenenergies in |2) but also the 
distribution of the coupling matrix elements A^. This distribution can be influenced by 
the degree of coherence in the wave packet \F) = F\0), initially formed by an external 
field. The motion in the isoscalar case is more coherent than in the isovector case 
since protons and neutrons move in phase. Therefore, the isoscalar state is expected 
to be more resistive against decay than its isovector counterpart ||. This tendency 
is also seen from Fig. 3 which displays the number of matrix elements A F ^ of given 
magnitude, connecting the state \F) and the states |2). In the cases referred to as (1) 
and (2) above, additional selection rules lead to large degeneracies and consequently the 
corresponding distributions sharply peak at zero. When the full residual interaction is 
used in generating the vectors |2) (case (3)) these selection rules are removed and there 
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are no vanishing coupling matrix elements any more. In the isovector case the wings of 
distribution seem to be consistent with a Gaussian as one would expect for a random 
process. In the isoscalar case, however, Ap^ remains more localized around small values 
and a Gaussian fit to the wings is unsatisfactory. 

The analysis of the distribution of the mixing matrix elements is consistent with the 
results seen in Figs. 1 and 2; the isovector strength distribution is more affected by 
the complex background. Another and perhaps the most interesting effect is that the 
isovector strength is distributed much more uniformly, not only in excitation energy but 
also in magnitude. A bit of imagination may even suggest a certain kind of self-similarity 
regarding the clustering and the relative size of the transitions. The picture in Fig. Id 



becomes reminiscent of a self-organized system at its critical state |L2], [H| where the 
equilibrium balance reduces the dimensionality. This observation finds confirmation in 
more quantitative terms. Fig. 4 shows the total number N of transitions of magnitude 
smaller than a given threshold value S t h, as a function of S t h- For the isovector resonance 
in the chaotic case (Fig. Id) we find, except for the largest transitions, a scaling law 
of the form N ~ Sf h (a ~ 0.50) (indicated by the straight line fit in the upper part of 
Fig. 4) which indeed signals a reduction of dimensionality. The two nonchaotic cases 
(Figs, lb and c) display a more complicated behavior. It is interesting to notice that a 
similar scaling applies also to the isoscalar resonance (lower part of Fig. 4) even though 
this is not so obvious from Fig. 2d. This time, however, the scaling interval is somewhat 
shorter and a ~ 0.48. 

The strength distribution can be considered as an attractor for the decay process 
starting out of equilibrium. It is nothing but the Fourier transform of the time correlation 
function (F(0)\F(t)) which, in the form of an envelope |14]], describes the process of 
gradual convergence to such an attractor and, thus, resolves it. This sets the parallel 
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to a procedure [15j which resolves the self-similarity and the scale invariance in classical 
chaos. This analogy provides further arguments for interpreting the above scaling law 
as another manifestation of 'l/f'-type behavior So far, such a behavior has been 
identified mostly on the classical level |17| for a variety of observables and models. In 
the present case of an 'avalanche' of the decaying giant resonance it is reflecting the 
chaotic properties of a strictly quantum mechanical phenomenon. 



The influence of chaos on the strength distribution is believed [JTSj, to manifest 



itself in a Porter-Thomas distribution [gOj of transition strengths maximizing the en- 
tropy of the strength distribution [IS]. Converting the appropriate (differential) Porter- 
Thomas distribution to a cumulative allows a direct comparison with our calculated 
distributions. The results are given by the thin solid lines in Fig. 3. They show the 
same scaling behavior with a = 1/2 exactly as can be easily seen from the analytical 
form of the Porter-Thomas distribution. One thus finds an impressive consistence. The 
strength distributions in Figs. 1 and 2 represent the coexistence of collectivity and chaos 
and, therefore, the amplified abundance of large components as compared to a purely 
random process is natural. Being more coherent, for the isoscalar resonance the onset 
of the scaling regime is delayed. It is also interesting to see that even a fully random 
process, as represented by the Porter-Thomas distribution, develops a trace of such a 
delay. The possibility of a similar two-phase behavior has been identified PI] in random 
cascade models. 



Acknowledgement 

This work was supported in part by the Polish KBN grants 2 2409 9102 and 2 P302 
157 04 and by NSF grant PHY-89-21025. One of the authors (S. N.) would like to 

7 



express his thanks to the Alexander von Humboldt Foundation for a fellowship. 



8 



References 

[1] T.A. Brody, J. Flores, J.B. French, P.A. Mello, A. Pandey and S.S.M. Wong, Rev. 
Mod. Phys. 53, 385(1981); 

R.V. Haq, A. Pandey and O. Bohigas, Phys. Rev. Lett. 48, 1086(1982); 
O. Bohigas, M.J. Giannoni and C. Schmit, Phys. Rev. Lett. 52, 1(1984) 

[2] T. Guhr and H.A. Weidenmiiller, Ann. Phys. 193, 472(1989) 

[3] J.J. Verbaarschot, H.A. Wiedenmiiller and M.R. Zirnbauer, Phys. Rep. 129, 
367(1985) 

[4] T. Srokowski, J. Okolowicz, S. Drozdz and A. Budzanowski, Phys. Rev. Lett. 71, 
2867(1993) 

[5] V.G. Soloviev, Nucl. Phys. A554, 77(1993) 
[6] V.G. Zelevinsky, Nucl. Phys. A555, 109(1993) 

[7] S. Drozdz, S. Nishizaki, J. Speth and J. Wambach, Phys. Rev. C49, No. 2 (1994), 
to appear 

[8] G.F. Bertsch, PF. Bortignon and R.A. Broglia, Rev. Mod. Phys. 55, 287(1983); 
S. Drozdz, S. Nishizaki, J. Speth and J. Wambach, Phys. Rep. 197, 1(1990) 

[9] B. Schwesinger and J. Wambach, Nucl. Phys. A426, 253(1984) 

[10] A. van der Woude, in Electric and Magnetic Giant Resonances, p. 99, ed. J. Speth, 
(World Scientific, 1991) 

[11] M.V. Berry and M. Tabor, Proc. R. Soc. Lond. A356, 375(1977) 

[12] P. Bak, C. Tang and K. Wiesenfeld, Phys. Rev. Lett. 59, 381(1987); Phys. Rev. 
A38, 364(1988) 



9 



[13] L.P. Kadanoff, S.R. Nagel, L. Wu and S. Zhou, Phys. Rev. A39, 6524(1989) 

[14] E.J. Heller, J. Chem. Phys. 72, 1337(1980); Phys. Rev. A35, 1360(1987) 

[15] S. Drozdz, J. Okolowicz and T. Srokowski, Phys. Rev. E48, no. 6 (1993) 

[16] P. Dutta and P.M. Horn, Rev. Mod. Phys. 53, 497(1981) 

[17] J.M. Carlson and J.S. Langer, Phys. Rev. Lett. 62, 2632(1989); 

K.L. Babcock and R.M. Westervelt, Phys. Rev. Lett. 64, 2168(1990); 
M. Ploszajczak and A. Tucholski, Phys. Rev. Lett. 65, 1539(1990); 
H. Takayasu and H. Inaoka, Phys. Rev. Lett. 68, 966(1992); 
L. Oddershede, P. Dimon and J. Bohr, Phys. Rev. Lett. 71, 3107(1993); 
B. Drossel, S. Clar and F. Schwabl, Phys. Rev. Lett. 71, 3739(1993) 

[18] Y. Alhassid and R.D. Levine, Phys. Rev. Lett. 57, 2879(1986) 

[19] M. Matsuo, T. Dossing, E. Vigazzi and R.A. Broglia, Phys. Rev. Lett. 
2694(1993) 

[20] C.E. Porter and R.G. Thomas, Phys. Rev. 104, 483(1956) 
[21] A. Bialas and K. Zalewski, Phys. Lett. 238B, 413(1990) 



10 



Figure Captions 

Fig. 1 The isovector quadrupole strength distribution in 40 Ca: (a) no coupling to the 
2p-2h subspace (b) no residual interaction in 2p-2h subspace (c) including only 
particle-particle and hole-hole matrixelements in the diagonalization of the 2p-2h 
subspace (d) diagonalization of the full residual interaction in the 2p-2h subspace. 

Fig. 2 Same as Fig. 1 but for the isoscalar quadrupole strength distribution. 

Fig. 3 The distribution of coupling matrixelements of the isovector and isoscalar wave 
packets \F) to 2p-2h excitations: (a) with no residual interaction in the 2p-2h 
subspace, (b) including only pp- and hh matrixelements in this space, (c) diago- 
nalizing the residual interaction fully. The full lines denote a Gaussian fit to these 
distributions. 

Fig. 4 The total number N of transitions of given strength below a threshold value S t h 
as a function of Sth- The open triangles refer to the case (b), the open squares to 
the case (c) and thick dots to the case (d) of Figs. 1 and 2, respectively. The solid 
line indicates the best fit to the later case. The thin solid lines represents the same 
quantity determined from a Porter-Thomas distribution. 
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